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theories in the limit where all kinematic invariants are large compared to particle masses. We 
use this relation to derive the two-loop QED corrections to large-angle Bhabha scattering. 
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I. INTRODUCTION 

Bhabha scattering, e + e~ — > e + e~, is one of the basic processes at electron-positron colliders. 
It has sufficiently large cross-section to be employed as a reference process for collider luminosity 
measurements. To determine the luminosity, one takes the ratio of the number of electron-positron 
pairs observed in the detector and the cross-section for Bhabha scattering computed theoretically. 
It follows that the accuracy of the luminosity determination is fully correlated with the precision 
of the theoretical description of Bhabha scattering. 

There are two kinematic regimes relevant for Bhabha scattering. Small angle scattering, 9 < 
few degrees, was employed at LEP and SLC for the luminosity determination, while large angle 
Bhabha scattering, 9 > 10 degrees, is used at flavor factories (BABAR, BELLE, BEPC-BES, 
CLEO-C, DAPHNE, VEPP-2M) for that purpose 0. Large angle Bhabha scattering can also be 
used to determine the luminosity spectrum at the International Linear Collider (ILC). 

Because of the importance of Bhabha scattering, the theoretical description of this process 
is quite advanced. In particular, several Monte Carlo event generators have been developed to 
describe the Bhabha scattering process |3|; as a rule, these programs correctly reproduce the 
0(a) QED corrections and some parts of higher order corrections enhanced by large logarithms. 
Further progress in the theoretical description of Bhabha scattering requires the computation of 
0(a 2 ) next-to-next-to-leading (NNLO) QED corrections and the incorporation of those corrections 
into Monte Carlo event generators. The virtual two-loop corrections to e + e~ — ► e + e~ scattering 
amplitude were calculated in Ref. in the approximation where lepton masses were set to zero. 
Unfortunately, this result is insufficient for the description of Bhabha scattering at flavor factories 
and the ILC where isolated leptons rather than "QED jets" can be observed. Furthermore, all 
existing event generators work with massive electrons. To include the NNLO results into these 
codes, it is thus necessary to keep the electron mass also in fixed order calculations. 

For all practical applications, having the logarithmic electron mass dependence of the NNLO 
corrections is sufficient. Corrections suppressed by powers of the electron mass are negligible, even 
at the smallest measured scattering angles. When performing the calculation of the cross section 
with a nonzero electron mass, one can thus expand the relevant diagrams in powers of the electron 
mass. Even so, the evaluation of the two-loop corrections remains a formidable task. Results for 
some of the necessary loop integrals were presented in Refs. but so far only 

the part of the cross section which involves closed electron loops has been evaluated 11 , l3, 3j- 



It is possible to avoid the evaluation of the diagrams in the massive case. Instead, one can use 



the massless result of Ref. 14|] and restore the mass dependence in the logarithmic approximation. 



For corrections that do not involve closed fermion loops, to which we will refer to as "photonic" , this 



was shown in Ref. 15[] where the phenomenologically relevant result for these 0(a 2 ) corrections 



to the large-angle Bhabha scattering was first reported. However, the approach of Ref. 15[ is 
somewhat complicated since it requires to transform infrared and collinear 1/e poles, inherent 
to massless amplitudes computed in dimensional regularization, to In A and lnm e terms in the 
massive amplitude where the photon mass A regularizes infrared divergences and the electron mass 
m e regularizes collinear divergences. 

In this paper we point out that a much simpler procedure for deriving massive amplitudes 
from massless ones exists since the two amplitudes are related by multiplicative renormalization 
factors. These renormalization factors can be deduced from the knowledge of massive and massless 
electron Dirac form- factors. The massive amplitude constructed along these lines has its infrared 
divergences regularized dimensionally and collinear divergences regularized by the electron mass. 
Beyond its simplicity, the advantage of this approach is that it can be directly applied to QCD 



whereas the method of Ref. 15j relies on the photon mass as infrared regulator. In a recent paper, 



Mitov and Moch have obtained a similar relation between massless and massive amplitudes [161 ]. 
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In fact, for photonic corrections our relation reduces to their result. In addition, our method allows 
us to also treat contributions involving massive fermion loops. 

We apply our method to compute the NNLO QED corrections to large-angle Bhabha scattering; 
we include both photonic corrections and contributions from closed lepton loops. This calculation 
is a nontrivial application of our method. In addition, it provides an independent check of the 
computations of Refs. U, l2, 13, 15| with which we find complete agreement 1 . We also derive the 
NNLO contribution from loops with leptons heavier than the electron, i.e. muons and tau leptons, 
which was not available in the literature. 

The paper is organized as follows. In the next Section we present our notation and discuss 
the perturbative expansion of the large-angle Bhabha scattering cross-section. In Section IIIII we 
explain the factorization formula considering the electron Dirac form factor as an example. In 
Section ITVl we apply the factorization formula to compute the NNLO QED corrections to Bhabha 
scattering. We conclude in Section |V] Some useful formulas are collected in the Appendix. 



II. NOTATION 

Consider the process e + {p\) + e~(p2) — > e + (ps) + e~(p4) for energies and scattering angles such 
that the absolute values of all kinematic invariants (jp\ +P2) 2 = s, (pi — p^) 2 = t and (pi — pi) 2 = u 
are much larger than the electron mass squared, s, \t\, \u\ 3> m 2 . 

We compute the Bhabha scattering cross-section perturbatively in the on-shell scheme where 
the fine structure constant a is defined through the photon propagator at zero momentum transfer. 
Neglecting corrections suppressed by powers of the electron mass, we write 



da 



a 2 ( 1 — x + x 2 



1 + 



a\ fa 
- Ml + - I °2 

IT J \TT 



5 2 + 0(a 3 ) 



(1) 



where x = (1 — cos#)/2 and 9 is the scattering angle. We take into account electron and muon 
contributions to photon vacuum polarization diagrams in the approximation s ^> m 2 ^ ^> m 2 . For 
phenomenological applications, the contribution of r leptons may be required; it can be obtained 
from the formulas below with the obvious modification — > m T provided that the high-energy 
condition s 3> m 2 is valid. 

The higher order corrections to the Bhabha scattering cross-section depend logarithmically on 
the mass of the electron. Also, in order to arrive at a physically meaningful result, we need to 
allow for soft radiation with the energy of each emitted photon below some value u; cut -C m. 
The perturbative corrections depend logarithmically on cj cu t- The corrections sensitive to soft and 
collinear physics are numerically enhanced relative to other corrections; it is therefore customary 
to separate out those corrections when presenting results for perturbative coefficients. The 0{a) 
correction in Eq. (pQ) is well-known [l8| ; in the limit of small electron mass it can be written as 



Sx = [ 4L soft + 3 + -N f 



hi 



111; 



(0) 



(2) 



where L so f t = ln(2o; cu t / V*) • We have introduced a label Nf to distinguish corrections due to closed 
lepton loops. In particular, including electron and muon loops corresponds to Nf = Ne + N^ where 
Ni is the number of leptons of the ith flavor. To obtain numerical results one has to set = 1 
and N,. = 1. 



1 Note that Appendix B of Ref. [13 ] contains a misprint fv^ . 
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The part of the one-loop correction that is not enhanced by the logarithm of the electron mass 
to center-of-mass energy ratio reads 

,2 



c(0) 



-4 + 4 In 



1 — x 



L 



2iV m£ 2vr^ 

soft In — | - 4 

6 mi 



-2Li 2 (x) + 2Li 2 (l - x) - + f(x) 



(3) 



The function /(x) is defined as 

f(x) = (l-x + x 2 y 2 | Q - |a; + |jx 2 - H x 3 + 1^4^ ^2 + ^ 3 _ 4x + 9^ _ 3 A 



+ 



3 x l 
-x 



-x 3 + x 4 ln 2 (x) + 



x 3 + ( 2 - 4x + -x 2 - x 3 ) ln(x) 



ln(l — x) 



7 



-1 + -x x z + -x 6 



x 4 ^j ln 2 (l - x) + Nf (- - |Vl - x + x 2 ) ln(x)| 



(4) 



The second term in the perturbative expansion of the Bhabha scattering cross-section is en- 
hanced by up to three powers of the logarithm of the electron mass. We write the NNLO correction 
as 



N 



An 3 



in- 



+ 4 2) ln 2 (4 
\mi 



+ 4 1} In . ., 
1 mi 



where 



c(2) 



8Lf oft +(12 + -iV / )L soft + -+iV / 



°2 



+N f 



•16 + 16 In 



8 , x 
- in 

3 1 - x 



3 

x 



(-1 



In- 



55 



V 3 l-x 18 



+ 4°\ 



Ni 



+ -11+1^ ln^ + ^; 



mi 



1 — x 
64\ 
~9~J 



-^soft + 



-28- V + 12 In ' 



+ 4/(x) 



8iV„ , m 



+6Li 2 (l-x) + 3/(x)+iV A( 



o m ,i 

- In 2 + - In 



In- 



mi 



L. 



soft 



mi 



1 — x 
93 5tt 2 



8Li 2 (x) + 8Li 2 (l-x) 
+ 6C(3) -6Li 2 (x) 



(5) 

(6) 
(7) 



m, 



281 



Nf ( -Li 2 (x) + — ) + N f g(x) + Nf [ -- + 



10 



l-x 9 

(2-x) 
9 ' 3(1 -x + x 2 ) 

and the function g(x) reads 

g(x) = (1 - x + x 2 ) -2 { (-x 4 - -x 2 — —x 3 + —x - - ) ln 2 (x) 
yv ; v ; IV3 4 12 12 3/ u 



ln(x* 



2 mf, 
3^ ln 4' 



+ ((1 - 2x) Qx 3 - X -x 2 + |) ln(l - x) + 
3 



37 56x 47x 2 67x 3 10x 4 * 



+ 



2 7 2 

-x 2 x + - ) (x 2 - x + 1) ln 2 (l - x) + 



6 

31 

18" 



10 2 31 , 
"^"X + —x + —x 



18 
31 
18' 



+ 



ln(x) 



£ - ^ ln(l 
9 9 J 



11 



+ — x H — x 
1 12 9 



1 



H — x 



23 
18 : 



-X 7T 



(8) 



Except for the muon contributions, the terms 6>p and S^ 1 ' were computed in Refs. 19j, l2Cll . |21 | 
and the term 4 in Eq. ([5]) was computed in Refs. 11 , 12, 13, 15]. We present the result of our 
computation of the term 4°^ below. 
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III. MASS FACTORIZATION 



We begin with the description of the method that we use to compute the NNLO QED corrections 
to the Bhabha scattering cross-section. The key to our approach is a factorization formula that 
relates massless and massive amplitudes for a given process; such a relationship between massless 
and massive amplitudes can be ejected because of tne well-known fact thai in physical gauges 
collinear divergences factorize into the wave function renormalization constants. 

We will explain the factorization theorem using the Dirac form factor and then apply it to 
Bhabha scattering but we stress that the same relation applies to arbitrary scattering amplitudes 
in QED and QCD in the limit where all particle masses are much smaller than typical momentum 



transfers. The Soft-Collinear Effective Theory (SCET) 231. \24. 1251] is an appropriate framework to 
analyze factorization properties of processes in this limit. This effective field theory is constructed 
by studying the perturbative expansion in QED or QCD and identifying those momentum regions 
in loop integrals that lead to singularities once the expansion of diagrams in small kinematic 
variables or masses is performed. These momentum modes are described by effective theory fields 
while the remaining contributions are integrated out and absorbed into the Wilson coefficients 
of the effective theory operators. The singularities relevant to our case arise when particles are 
soft or have momenta collinear to the external momenta. The effective theory description of the 
electric current requires two different collinear fields which interact via soft exchanges. To explain 
the structure of the result, we first consider the electron Dirac form factor in the limit where the 
electron momenta fulfill Q 2 = — (pi — P2) 2 3> p\ ~ p\ ~ m 2 . At leading power in the effective 
theory, the vector current = ipjuip takes the form 

= J dsdtC v (s,t) [&W 2 ](sn 2 )^ [Wtii](tn x ) + OiQ- 1 ). (9) 

Here, ni and ri2 are the light-like reference vectors in the directions of p± and P2, respectively. The 
conjugate vectors fii and n<i point in the opposite directions and fulfill ni • ni = T12 ■ n<i = 2. These 
reference vectors must be chosen such that fii ■ pifi2 ■ P2 = Q 2 + 0(m 2 ). The collinear electron 
fields £1 and £2 are multiplied by light-like collinear Wilson lines 

t ? \ 

Wi(x) = exp lie dsrii ■ A Cj i(x + sni) . (10) 

V -L I 

The Fourier transform of the Wilson coefficient 

C v (Q 2 ) = C v (ni ■ pin 2 -P2) = J dsdt C v (s, t ) e isfl ^' m ^ (11) 

depends only on the hard scale Q 2 , but is independent of the electron mass. To obtain Cv(Q 2 ), 
we perform a matching calculation. The simplest way to do the matching is to use dimensional 
regularization and to calculate the on-shell form factor in the massless theory. In this case all 
loop diagrams in effective theory vanish and the bare Wilson coefficient Cy(Q 2 ) equals the on-shell 
Dirac form factor F(Q 2 ) of a massless electron. The massless on-shell form factor has infrared 
divergences, which show up as poles in 4 — d = 2e. These poles correspond to ultra-violet diver- 
gences in the effective theory. Since the Wilson coefficients are independent of the small electron 
mass, the difference between massive and massless amplitudes can only arise from matrix elements 
of operators in the effective field theory. Off-shell Green's functions in the effective theory get 
contributions from soft interactions between external legs and collinear interactions in each sector. 
At leading power, soft photons have eikonal interactions with collinear fields; only the • A s com- 
ponent of the soft photon field interacts with collinear electrons moving in the ith direction. These 
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interactions can be removed by field redefinitions [2J 

&(s) = 5 4 (x„)ef ) (x), A?(x) = Si(x-)4 0)fl sl{x-), (12) 
where x_ = ^(n, • x)n^ and the soft Wilson line reads 

I ? \ 

Si(x) = exp \ie dsrii ■ A c ^(x + sni) . (13) 

V -co / 

The current operator takes the form 

V, = fdsdtC v (s,t) [ef Vf ]( S n 2 ) 7 ^|(0)5 1 (0)[^ 0)t ^ ) ](tn 1 ) +0(Q- 1 ). (14) 



After the field redefinition Eq. (|12p . there is no interaction between the different sectors of the 
theory. The matrix elements of the current operator factorize into collinear matrix elements for 
each direction, called jet-functions, and a soft function, which is given by the matrix element of 
the soft Wilson lines. This factorization of Green's functions into hard-, jet- and soft functions at 



large momentum transfers is a well known property of gauge theories [26j, l27|, |28| . 

For on-shell matrix elements the situation is especially simple because in the massless case 
jet and soft functions are trivial since effective theory loop diagrams are scaleless. For massive 
electrons, the collinear matrix elements are functions of the lepton masses, while the soft function 
also depends on the hard momentum Q. We therefore write 

F(Q 2 , {m 2 }) = Zj({m 2 })S(Q 2 , {m 2 })F(Q 2 ) + 0(m 2 /Q 2 ) . (15) 

The relation between the massive form factor F and the massless form factor F simplifies when 
only photonic corrections are considered. In that case higher order corrections to the soft matrix 
element in the massive theory are given by scaleless integrals and therefore vanish. This can be seen 
diagrammatically before performing the somewhat formal decoupling transformation in Eq. (|12p . 
Hence, we conclude that for photonic corrections the soft function in Eq. (|15p is equal to one to 
all orders in QED perturbation theory and can be dropped from the right hand side of Eq. fjl5jl . 
The relation between massless and massive form factors obtained in this way coincides with the 
relation discussed recently in Ref. (ill . It is also consistent with the one-loop relation obtained in 



Ref. |30|]. 



However, a non-trivial soft function appears once vacuum polarization diagrams with massive 
particles are considered. In that case, it is easy to see that the soft momentum contribution in 
the massive theory is not a scaleless integral and therefore does not vanish. Moreover, it exhibits 
non-trivial dependence on the hard scale Q. We write the soft matrix element S(Q 2 ,{m 2 }) as 



S(Q 2 ,{m 2 }) = l+ ]T 5S(Q 2 ,m 2 ), (16) 

where 



i=e,fi 



W?,*S,*) = J (g, ¥i . ^ k) p .n(^). (17) 

In Eq. (|17p oq stands for the bare QED coupling constant. The vacuum polarization function 
U(k 2 ) in Eq. (jUJ) is defined as 



,2 2, _ (-i) r d d i 



*n(/c ,m i ) = jj—[jk 2 ~ J (2^p Tr 



1 — mi f + If, — mi 



(18) 
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The easiest way to compute the soft matrix element is to employ a dispersive representation for the 
vacuum polarization function IL(k 2 ). Evaluating the integral in Eq. (117p in d = 4 — 2e dimensions, 
we obtain 



5S(Q 



2 2 

mi, 



Ni) = JViagni: 




1 t 5 
12? 36e 



7 



7T 



27 72 



(19) 



where ao = ao/ 71 " e _7e (4-7r) e and 7 is the Euler constant. Eq. (|15|) provides a relation between 
closed lepton loop contributions to massive and massless form factors and allows a derivation of 
A^e iAt -dependent 0(a 2 ) corrections to Bhabha scattering from the massless results of Ref. [141 ] . 

To determine the square of the jet-function Zj to NNLO in QED, we use Eq. (|15p and divide 
the ratio of dimensionally regularized massive and massless form factors by the soft matrix element 
Eq. (|17p . The expansion of the massive vector form factor in the limit Q 2 S> m 2 throu gh Q (a 2 ) 



can be found in Refs. [3l|, ]32j, ]33j . The massless result can be obtained from Refs. [34j, ]35j]. We 



express the jet function through the bare QED coupling constant an- We find 



Zj = 1 + a m, 



-2c 



1 1 IT 2 7T 2 

2? + 47 + 24 + 1 + e 2+ 48 



-a\m e 4£ 



1 1 

+ -s 



8e 4 



8 24 



+ 



1 





209 5vr 2 " 
432 + 144 



561 61vr 2 
128 + 192 




_C(3) vr__ 
12 320 

vr 2 2C(3) 




ic(3)-^ln(2) 



3379 



197T 2 C(3) 



2592 216 
7^ 2 ((3) 



1296 + 54 




l ln 3^ 
36 



mf 



25 , 

H In 

72 



2^ 
mi 



+ 




193 7T 

216 + ' ~ 



s 2 
ae , a e 



(20) 



Setting iVj = N„ = 0, we find the agreement with the result of Ref. la]. The independence of 
the jet function Zj of the hard scale Q 2 is an explicit demonstration of the mass factorization to 
two-loop order. 

Before turning to Bhabha scattering, we want to address a subtlety concerning Eq. (]15p . This 
relation relies on the assumption that only hard, collinear and soft momentum modes are relevant 
in the effective theory computation. However, as was explicitly shown in Ref. [29j], this assumption 
is invalid for some diagrams that contribute to the form factor. A particular example discussed in 
Ref. 29J] is the contribution of so-called ultra-collinear modes to a two-loop planar vertex diagram. 
Let us stress that the relevance of these modes for the full form factor computation would invalidate 
the factorization formula Eq. (I15p since this mode induces additional dependence on the hard scale 
Q that is not associated with the hard or soft region. It is therefore gratifying to observe that 
the ultra-collinear modes discussed in Ref. [29| are not relevant for the form factor since their 
contributions cancel out. For example, at the two-loop level, the ultra-collinear contribution to 
a planar vertex diagram cancels exactly against a similar contribution to the non-planar vertex 
diagram making the full form factor independent of it. 



2 We need the 0(e 2 ) terms in the one-loop contribution to the massive form factor which are not provided in Ref. [3lj . 
However, it is simple to compute these terms and the corresponding result is given in Appendix A. Note also that 
with the normalization adopted in Eq. (24) of Ref. [3l| the MS result for the contribution of massive quark vacuum 
polarization to the heavy quark form factor coincides with the result in the on-shell scheme. 
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IV. BHABHA SCATTERING 



We can now apply the factorization formula, established for the electron Dirac form factor in the 
previous Section, to Bhabha scattering. To get the scattering amplitude M in which the electron 
mass is used as a regulator of the collinear singularities, we only need to multiply the massless 

1/2 

amplitude M by the square root of the jet function Zj for each electron and positron leg and by 
the product of soft functions that account for soft exchanges in s,t and u channels. To obtain the 
soft matrix element, we may use Eq. ()19j) . While Eq. (|19p is relevant for the space-like electron 
form factor, we need to generalize it to describe soft exchanges in the s- and u-channel. We write 

M({ Pi },{m 2 }) = Z 2 ({m 2 })M({ Pl })S(s,t,u) + 0(m 2 /Q 2 ), (21) 

where the soft function S(s, t, u) is given by 

S(s, t, u) = I 1 + 2 J2 E 8S(Q 2 ,ml XguNi) J (22) 

and the sum goes over Q 2 = —s, —t, —u with A s = A^ = 1 and X u = —1. The change Ni — > — N{ 
shown in Eq. (|22p accounts for the required change in the overall sign in 5S in the it-channel. 
We can now use the two-loop result for the Bhabha scattering amplitude computed in the 



massless limit 14] and employ Eqs. (|20p . (|2ip and (|22p to obtain the scattering amplitude in which 
collinear divergences are regularized by the electron mass and infrared divergences are regularized 
dimensionally. With the massive scattering amplitude at hand, the computation of the two-loop 
QED corrections to the large-angle Bhabha scattering cross-section becomes straightforward; what 
we need in addition is the cross-sections for the inelastic processes e + e~ —* e + e~ + 717, with n = 1 
and n = 2, in the soft photon approximation. We write the perturbative expansion of the Bhabha 
scattering cross-section as 

^ = exp (-F soft ) Z) \S\ 2 - (da + -Aa\ + (-Y (d< xl + da v 2 ) + 0(a 3 )) . (23) 

In this formula, Zj is the square of the jet-function given in Eq. (|20p ; -F so f t describes soft photon 
radiation which, in case of QED, is known to factorize and exponentiate. Up to an overall factor 
a 2 / s the quantity do"o is the tree level cross section in d-dimensions and daf denotes the one- 
loop virtual contributions. At two loops, there are two types of virtual corrections: the quantity 
da?? contains the interference of the two-loop amplitude with the tree level amplitude, while da\ xl 
describes the interference of one-loop amplitude with itself. These contributions are to be computed 
with massless leptons. The massive result for the virtual corrections is then obtained by multipl yin g 
the massless result with Zj\S 2 \. The massless cross-sections and da% can be found in Ref. [14] , 



while da\ x i can be obtained from Ref. 36f|, as we explain below. Note also that results of Refs. [14l . 
361 ] are written through the QED coupling constant a renormalized in the MS scheme. It is for 
this reason that a appears in Eq. (I23p . 

In Ref. [i^] the interference of the one-loop amplitude with itself was obtained for quark-quark 
scattering in QCD. To extract the QED piece relevant for Bhabha scattering from these results, 
we need to analyze the color algebra; such an analysis shows that da\ xl can be obtained from the 
computation of Ref. [36] by taking the — > limit of the QCD result, where A^ is the number 
of colors, and subtracting from it suitably weighted products of the one-loop Bhabha scattering 
cross-section da\ and the electron Dirac form factor in the massless approximation. Note that 
the divergent terms in daf xl can be obtained from Catani's decomposition [37| of the one-loop 



9 



scattering amplitude for e + e~ — > e + e~ and the O(a) correction to the Bhabha scattering cross- 
section in dimensional regularization derived in Ref. [141 ] . For this reason we only present the finite 
part of daf xl ; it is given in the Appendix B. 

Finally, the description of soft radiation in Eq. ()23f) is encapsulated in the function i^ so ft; we 
require this function for massive electron-positron scattering. Since i^oft describes the emission of 
a real photon, it is simplest to evaluate it in the on-shell scheme. In doing so, the effects of vacuum 
polarization contributions are absorbed into the coupling constant and do not need to be evaluated 
explicitly. The function F so f t is determined by the integral 



F, 



soft 



-Air 1 



fco<Wcut 



ddk J J, 

(2n)^2k ' 



(24) 



where the soft current is given by 



Pi-k' 



(25) 



qi is the charge of the particle i in units of the positron charge and A = ±1 for the incom- 
ing(outgoing) particle, respectively. Writing 



7 W = V Pi ' Pj 

t£ (Pi ■ k)( Pj ■ k) 



qiqjXiXj + 



ml 



(pi • k) 2 ' 



(26) 



we observe that two types of integrals are required for the evaluation of F so f t . We give expressions 
for these integrals below neglecting all the terms that are suppressed by powers of the electron 
mass. 

The first integral depends on the relative momenta of two charged particles. We find 



d d-l k 



Pi ■ Pj 



'l 



{2TT) d - 1 2k { Pi ■ k)(pj • k) 



-K 



(2w. 



cut / 



-2c 



2e 



In 



1 



1 + Xi 



7T 



+ e o^m + Li 2 (-Xij) + — + € 



r3 
6 



L m + 2C(3) +Li 3 (-xij) 



(27) 



where = (1 + cos^)/(l — cos 9ij), 9ij is the relative angle between the three-momenta of the 
particles i and j and the normalization factor A4 reads 

T(l-e) 



47r 2 (47r)- e r(l - 2e)' 



(28) 



Also, we introduced L m = ln(mg/s) to denote the collinear logarithm. The second integral required 
to describe the soft radiation reads 



d^k 



ml 



fco<OJ cu t 



(2vr) d - 1 2fc { Pl -kf 



-K 



(2w. 



cut / 



-2c 



2e 



l-eln^ + -ln z ^ 



(29) 



where we used the on-shell condition pf = m\. 

We now have everything in place to calculate the NNLO QED corrections to Bhabha scattering. 
We substitute all the necessary ingredients into Eq. (|23|) . To combine the different pieces, it is 
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simplest to first express the on-shell coupling constant appearing in the soft radiation exponential 
in Eq. (|23l) through the MS coupling constant. In d-dimensions, the relevant relation reads 



a = a(/i) /j, 2e I 1 



2aQ) 
3vr 



E 

/=e,M 



2 M , ^ 



m/ 24 



(30) 



After adding the real and virtual parts, all infrared divergences cancel and we can use the four- 
dimensional relation between the MS and the on-shell fine structure constants 



a(fi) = a < 1 + ( — 



it J V 3 /// 



e 



+ - 



16 
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to obtain the expansion of the Bhabha scattering cross section through NNLO in terms of the on- 
shell QED coupling constant. Upon doing so, we reproduce the formulas for radiative corrections to 
Bhabha scattering shown in Section [II] and obtain an explicit expression for 4 , which is presented 
below. We write 82 in the following way 
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where we separate photonic corrections and the corrections caused by closed lepton loops. We now 
present the results for these terms separately. For the photonic corrections we find 
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The contribution of diagrams with a single electron or muon vacuum polarization insertion is 
described by S^' 2 ^ , S^' 3 ^ • We obtain 
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The corrections with two insertions of the closed lepton loop read 
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Because the above expressions are rather lengthy, we have included them in electronic form in our 
submission to the arXiv. 



V. CONCLUSION 

We presented a novel relation between massive and massless scattering amplitudes in QED 
valid in the limit when all kinematic invariants are large compared to masses of particles that 
participate in the scattering process; for quenched QED, our result agrees with a similar relation 
between massive and massless scattering amplitudes discussed recently in [3] . We used this relation 
to derive the NNLO QED corrections for Bhabha scattering confirming earlier results on photonic 
and electron loop corrections of Ref. U, E3, 13, We also obtained NNLO contributions to 



Bhabha scattering due to muon and tau vacuum polarization loops that were not available in the 
literature. 



In variance with the approach discussed in Ref. [15(, our method is directly applicable to QCD. 
A potentially interesting application is the computation of the NNLO QCD virtual corrections to 
heavy (e.g. b) quark production at moderate to large momentum transfers at the Tevatron and 



13 



the LHC using the two-loop matrix elements for the gg — > qq process computed in Ref. 38)] for 
massless quarks. 
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Note added While we were finalizing our paper, Ref. [39J appeared in which the muon-loop 
contribution to Bhabha scattering was evaluated. After the authors of Ref. [39] corrected one of 
the master integrals their result agrees with ours. 
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APPENDIX A: THE ONE-LOOP MASSIVE DIRAC FORM FACTOR 



The one-loop massive Dirac space-like form factor, in the limit Q 2 3> m 2 , reads 
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where L = ln(Q 2 /m 2 ) and a(m) is the MS fine structure constant evaluated at the scale m. 
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APPENDIX B: NNLO RESULTS 

Here, we give the finite part of the 0(a 2 ) correction to the cross-section due to the interference 
of the photonic one-loop amplitude with itself. The divergent part can be obtained using Catani's 



results [37j] 
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